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'^ I Motivated by the indication that both the atmospheric and the solar neutrino 

y—^ ■ puzzles may simultaneously be solved by (vacuum as well as matter-induced 

'"^ ' resonant) oscillations of two generations of neutrinos with large mixing, we 

Qs . have analyzed the data on the atmospheric and solar neutrinos assuming 

"^ \ that all three neutrinos are maximally mixed. It is shown that the values 

CIhI of A?7i^ obtained from the two-generation analyses are still valid even in the 

o ' three-generation scheme, i.e. the two puzzles can be solved simultaneously if 

^ ' ^^^31 — 10~^ eV^ for the atmospheric neutrinos and A?n,2i — 10~^° eV^ for 

• • . solar neutrinos in the maximally mixed three-generation scheme. 
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I. INTRODUCTION 

One of the several outstanding issues that we have faced so long in particle physics is that 
of neutrino masses and mixing (see, for example, Refs. [|l],3). Studies of solar neutrinos have 
been very promising for some time and already we have preliminary indications for positive 
answers 0,^,|5|J^. In this case, however, one has to rely on the solar model 0, which is as 
good as it can be but it may not be, as yet, perfect. For this reason, many suggestions have 
been proposed on solar model-independent ways of analyzing solar neutrino data when more 
data will become available in the future P,^. 

Another prominent source of neutrinos that has been investigated in the past is atmo- 
spheric neutrinos [|l^,|ll|,|l^,|l^,|l^,|l^,|l^ . In this case, although the absolute fluxes of u^ (i/^) 



and z/g (ue) are uncertain up to 30% [|r^, if one takes the ratio of neutrino-induced muon and 
electron events both for data and MC simulations (in fact the ratio of the ratios), the results 
are almost model independent with errors of about 5% |]T^ . In this sense, the atmospheric 
neutrinos may be better suited for extracting information on the neutrino parameters such 
as mass and mixing angle. 

Recently the Kamiokande collaboration |T^ has analyzed more than 200 events with 
energies in the multi-GeV range and updated the previous sub-GeV data. Based on these 
data they have concluded that both data can be explained by either z/g-t'^ or z/^-z/^ oscilla- 
tions. However, they have used explicitly the two generation formulations. The results of 
their analyses under the assumption that the anomalies are due to neutrino oscillations are 
summarized as being consistent with sin^(2'(9) ~ 1 and Am^ ^ 10~^, where f? and Am^ are, 
respectively, the I's-^/i mixing angle and m"^ — ml, and the z^^-z^t mixing angle and m| — m2 
for the z/e ^ z/^ and u^ ^ z/^ interpretations. (Whenever we mention p^, v^ and Vr we also 
infer the cases of z^e, z/^ and z/,-.) 

We also recall that in the case of solar neutrinos there are basically two popular ways to 
explain the so-called anomaly (under the assumption that the SSM is correct) with neutrino 



oscillations, one with vacuum oscillations and the other with the MSW effect [113]. Both 



cases contain solutions with almost maximal mixing of the neutrinos involved |T9| , ^ . It 
is to be stressed here that for reaching the respective conclusions in all these analyses two 
generations of neutrinos have been used. 

Now we may ask an obvious and logical question: Do these results remain valid even 
when the third generation neutrino is included? After all, we do know that there exist 
three neutrinos associated with e, /i and r. However, a three-generation analysis of any 
phenomenon is known to be very complicated and tedious, and often it is hard to see the 
physics involved [ ^ . 

In this work, we analyze the atmospheric and solar neutrinos using the three generations 
of neutrinos that are maximally mixed to begin with. Although this is, admittedly, an 
extreme case, it is a very good approximation in the event of large neutrino mixings, and, 
moreover, the algebra and analysis of the data become very simple and the physics becomes 
transparent. Furthermore, this does not imply a loss of generality if mixings turn out to be 
reasonably large because the qualitative results still remain valid. 

The three-generation neutrino oscillation with maximal mixing has also been discussed 
in the past [^] without taking into account the energy and distance dependence of the 



oscillating terms. In this case the survival probability for the solar neutrinos, P{i'e -^ ^e)^ is 



1/3 at all energies. In Ref. [^ also the atmospheric neutrino anomaly has been addressed 
by assuming that all Am^- are in the range (0.5 — 1.2) x 10~^eV^. 

Our main purpose in this paper is to see if and how the two generation results are 
modified when the third neutrino generation is introduced by taking into account the energy 
and distance dependence of the oscillating terms. 



II. MAXIMAL MIXING AMONG THREE NEUTRINOS 



We consider three neutrinos to be maximally mixed if their survival probabilities are 
all equal and for distances much larger than the oscillation lengths all the three survival 
probabilities average to 1/3. This can be accomplished with the following mixing matrix 
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The presence in U of the cube roots of unity 1, e^'^*/^, e^'^*/^ is necessary for U to be unitary. 
It follows that in the maximal mixing case CP must be violated. Using the mixing matrix 
!TD we obtain the survival probabilities 
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Here L, E and Attz?- 
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are in units of meter, MeV and eV^, respectively. On 
the other hand, the probabilities of transition between different neutrino fiavors are not 
identical. Instead, we have 
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It is interesting to notice that the mixing matrix for maximal mixing is not unique. 
In fact, all the matrices which are obtained from f/ by multiplying any number of rows or 
columns by a phase give the same oscillation probabilities 



), (2]^) and ( p. 51) . Furthermore, 
V'^ gives the same survival probabilities (|2.2| ), but all the phases 7r/3 in the transition 
probabilities ( |2.4| ) and ( p.5| ) change sign. Since the definition of maximal mixing requires 
the survival probabilities (|2.2| ), it follows that both V and [/^ give maximal mixing, but they 
produce two different phenomenologies for the transition probabilities. 



III. ATMOSPHERIC NEUTRINOS 

As in many previous works, our study of atmospheric neutrinos involves the ratio 
(z/^ + i^ij.) I [ye + ^e)- In actual experimental situations, this ratio manifests as an observed 
ratio of muon-like and electron-like events. Although the issue of whether or not the ob- 
served ratio is indeed small compared with the Monte Carlo (MC) simulated ratio is still 
not completely settled |lT0| , p!Tl , p!2| , p!3| , p!^ , p!5| , p!6| , the Kamiokande group recently presented two 
sets of impressive data, one involving the electron and muon events in the sub-GeV energy 
range and the other in the multi-GeV energy range. The Kamiokande collaboration found 



0.60^Q05 ± 0.05 for the sub-GeV energy range , 



R (/^/e)Mc 1 0.57l|];g^ ± 0.07 for the multi-GeV energy range 



where R (/^/e) and R (Ai/e)j^(-. are the ratios of yU-like and e-like events observed and 
calculated, respectively. In addition, they studied the dependence of the above ratios on the 
zenith-angle Q in the range — 1 < cos^ < 1, divided in 5 bins. The zenith- angle dependence 



of R'^^P for the sub-GeV data is rather flat (see Fig.6 of Ref. ||T2[)) whereas that of the multi- 
GeV data appears to deviate from being fiat (see Fig.4 of Ref. [1^1). Our analysis will be 
based only on the Kamiokande data. 

Let us consider the zenith-angle dependence of the Kamiokande data. For each bin i 
(i = 1, . . . , 5 for (cos 6*)^^^ 5 = —0.8 , —0.4 , 0.0 , 0.4 , 0.8, which correspond to the averaged 
distances (X)i=\ 5 = 10) 230 , 5, 157 , 852 , 54 , 26 Km), the number iVe, of e-like events, the 
number N^i of /i-like events and the ratio of ratios Rj are given by 

iV,, = ^ iV,^^P,(z/e - Ve) + PT V^iV,"^^ P^{y , - ^e) , (3.2) 

vPr 



^,^ = -^ N^"" P.(^e - ^,) + PT V^iV2^P.(^. - ^.) , (3.3) 
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p _ N,jNf^^ _ P.{^, ^ .,) + Pi^ [N^^'lN^^y PAve -. V,) 
''-Ne./N^- '^ P.(..^..)+pK(iV--/iVMC)p^(,^^,^) • ^'-^ 

Here iV^^ and iV^^ are, respectively, the number of e-like and /i-like events in the bin i 
predicted by MC simulations, without neutrino oscillations, after passing through the same 
analysis chain as the data. The factors px and pr are the normalization coefficients that 
take into account the systematic error of the MC calculation of the total number of e-like 
and p-like events (30%) and the systematic error of the MC calculation of the p/e ratio 
Vo for the sub-GeV data and 12% for the multi-GeV data), respectively. We inserted 



these two factors in Eqs. (|3.2| ) and ( |3.3| ) in such a way that the ratio of ratios Rj does not 
depend on px, i.e., only the coefficient p^ appears in the expression (|3.4| ) for the ratio of 
ratios as an overall multiplicative factor and the ratio [Nj^'^ /N^'^j is always multiplied 
by pr. As a consequence, the error of R, does not depend on the systematic error of the 
MC calculation of the total number of e-like and p-like events, but only on the systematic 
error of the MC calculation of the p/e ratio, which is much smaller. In Eqs.(|3.2|)-(|3.4|) we 



neglected the fact that not all e-like and /i-like events are produced by u^ ipe) and z/^ (z/^), 
respectively. However, the purity of the v^ and v^ contributions to the selected e-like and 
/i-like events is estimated by the Kaniiokande Collaboration to be higher than 90%. We wish 
to emphasize here that our three generation formalism takes into account simultaneously 
the two neutrino oscillation channels v^ ^ z/g and z/^ — > z/,-. This is not the case for two 
generation analyses, where the two oscillation channels are considered separately, leading to 
two possible interpretations and results, as done in Ref. |T2|. 



Since we have L < 13 x 10^ Km, E > 0.2 GeV and Amli ~ 10^^° eV^ (as will be seen 
later from the solution of the solar neutrino problem), ^21 is well approximated by zero and, 
because of the assumed hierarchy ms > 777-2 > mi, we expect /C32 — /C31, so that the relevant 
transitions probabilities are 

4 

P{l^e -^ Z^m) = Pil^fi -^ ^e) = P{.^t, -^ J^t) = P{.J^e -^ ^r) = - siu^ h2 ■ (3.5) 

y 

We obtained the values of N^'~' and N^*^ from the MC data for the e-like and /x-like 
events presented in Fig. 2 of Ref. []TU[ for the sub-GeV data and the MC data presented in 



Fig. 3 of Ref. [Q for the multi-GeV data. The quantities Pi(z^e — >■ i^n) and Pi(z/^ — >■ z/g) are 
the z/g —>■ v^ and u^ — *■ Ve transition probabilities averaged over the neutrino energy spectrum 
and the zenith-angle width of bin i. These probabilities depend on the value of Attt-Ij and 
they are different because initial z/g and z/^ have different energy spectra. These energy 



spectra are given in Fig.l of Ref. |]TT| for the sub-GeV data and the MC data presented in 
Fig.2 of Ref. H for the multi-GeV data. 

Before we present the results of our full analysis, we discuss the following simple but 
interesting possibility. As suggested by the zenith-angle distribution of the sub-GeV data 
and to some extent (within error bars) by the large-angle data points of the multi-GeV data 
with 2 = 1, ... ,4, corresponding to {cos 6) from -0.8 to 0.4, R'^^p is roughly consistent with 
being fiat. This may be interpreted as due to washing out of the distance dependence in 
the argument of sin^ {1.27LAml2/E) (the averaged distance (L) ranges from 10,230 Km to 
26 Km for the sub-GeV data and from 10,230 Km to 54 Km for the multi-GeV data). That 
is, the oscillations are in the rapid oscillation region with sin^ {1.27LAm22/E) ~ 1/2 and 
P(z/q, -^ Vjs) ^ 2/9 for a ^ (3. Under this assumption, Eq.( p.4|) reduces to 

R. = PK \ '" ^ ^' ' - ^ . (3.6) 

Since we have N^^/N^^ . ~ 1.57 for the sub-GeV data and N^^/N^^ 



2.98 , 2.29 , 2.08 , 2.39 for the multi-GeV data, and taking into account the systematic 
error of pr, we find Ri=i,...,5 = 0.83 ± 0.04 for the sub-GeV data and Rj=i,...,4 = 
0.62 ± 0.04 , 0.70 ± 0.05 , 0.74 ± 0.05 , 0.69 ± 0.04 for the muhi-GeV data. These values 
are to be compared with (R^""^) = 0.601o;q7 for the average of the sub-GeV data and 
R^^^P^ ^4 = 0.29l|J:J? , 0.46t[}:}^ , 0.5llg:J? , 0.63l|]:?^ for the multi-GeV data. It appears that 
the idea of rapid oscillations is not too inconsistent with the observations. 

Let us now discuss the results of a full analysis of the Kamiolcande data taking into 
account the oscillation effect. For different values of Am|2 we have fitted separately the 



sub-GeV data of the ratio of ratios presented in Fig. 6 of Ref. |T2[ and the multi-GeV data of 
the number of e-hke and fi-like events presented in Fig. 3 of Ref. 0]. For the fit of the sub- 
GeV data we used the theoretical formula for the ratio of ratios given in Eq.( p.4| ) and for the 
fit of the multi-GeV data we used the theoretical formulas for the number of e-like and /i-like 
events given in Eqs.( |3.2| ) and ( |3.3| ). For each value of Aml2 we have calculated the minimum 
value of the x^ with respect to the normalization coefficients px and pr, whose difference from 
unity contributes to the x^ with errors cr(pR) = 9% for the sub-GeV data and 0"(pt) = 30%, 
o'(Pr) = 12% for the multi-GeV data. Some results are depicted in Fig.|I| together with 
the experimental data. As can be seen in Fig.|T], the sub-GeV data and multi-GeV data 
can be best fitted (solid lines) with Aml2 = 1-3 x 10~^eV^ and Aml2 = 3.0 x 10"^ eV^, 
respectively. In Fig.| we have plotted the x^ oi the fits as functions of Am^g- As mentioned 
already, the best fits of the sub-GeV (FigJ^a) and multi-GeV (Fig.0b) data are obtained for 
Am|2 = 1.3 X 10"^ eV^ (x^ = 4.9 with 3 degrees of freedom) and Amfa = 3.0 x 10"^ eV^ 
(x^ = 8.1 with 7 degrees of freedom), respectively. The separate fits of the the sub-GeV and 
multi-GeV data are acceptable at 18% CL and 32% CL, respectively. The 90% CL allowed 
range of Aml2 is determined by x^ < Xmin + 2-7, that gives Am'^2 > 8.1 x 10"'^ eV^ from the 
sub-GeV data and 2.6 x lO'^eV^ < Am|2 < 6.3 x lO'^eV^ from muhi-GeV data. We see 
that these allowed ranges are perfectly compatible. Fig.^ gives the x^ of the combined fit 
of the sub-GeV and multi-GeV data. Here there is a minimum for Am'^2 = 2.0 x 10~^eV^, 
very close to the two independent best fits of the sub-GeV and multi-GeV data. Since 
Xmin = 13.8 and we have 11 degrees of freedom, the fit is acceptable at 24% CL. As can be 
seen from Fig.|^, the simultaneous fit of both the sub-GeV and multi-GeV data allows to 
constrain the value of Am^g to lie in the interval 8.9 x 10~^ eV^ < Am^g < 5.6 x 10~^ eV^ at 
90% CL. This result agrees well with those given in Ref. [12], obtained from two generation 
analyses. We wish to add here that the values of Am'^2 larger than 9.0 x 10~^ eV^ and larger 
than 2.0 x 10~^eV^ are excluded at 90% CL by the Kurchatov |2^ and Gosgen [^ reactor 



neutrino oscillation experiments, respectively. If combined, we are left with a very restricted 
allowed range for the value of Ar7i|2, which is about 10~^eV^. 

IV. SOLAR NEUTRINOS 

We now proceed to check if the three-generation maximally mixed scheme discussed so 
far is also capable of solving the solar neutrino problem. Since we assume Am|]^ ~ Am'^2 ^ 
10~^°eV^ as obtained in the atmospheric neutrino problem, we have, from Eqs.( p.4|) and 
T5D, 



2 9 9 

and, from Eq.( |2.2| ), 

5 4 



sin^A;2i, (4.1) 



P{u, ^u,) = ^-^ sin^ k2i . (4.2) 

Note that the maximal value of the survival probability of z/g is (5/9) = 0.56. This implies 
that, since the Standard Solar Model predicts the event rate for the Gallium detectors to be 
132lg SNU, the observed event rate must be 

6 



Sexp < 80 . (4.3) 

If the Eexp turns out to be greater than the above, this maximal mixing model (as an 
explanation with vacuum oscillation) is ruled out (or Am'^2 ^^^ ^'^ii niust take values 
such that sin^ ^32 and sin^ /C31 cannot be approximated by 1/2). The latest average of the 
GALLEX IQ and SAGE [Q event rates is Sexp = 74 ± 9 which is still consistent with Eq. 

The results of the calculations for the Kamiokande, Homestake and Gallium experiments 
are shown in Fig.^. These calculations are based on the standard procedure, with the 
neutrino fluxes predicted by the Bahcall and Pinsonneault SSM ^7^, the detector cross 
sections given in Refs. 0J2^, and the probabilities given in Eqs.( |4.1| ) and ( ^.2| ). 



In the figures, the horizontal solid lines are the central values of the experiments with 
the dashed lines representing one standard deviation. The oscillating solid lines represent 
the theoretical values, with the dotted lines indicating one standard deviation resulting from 
the errors in the SSM. We can see that in the neighborhood of Am'^i — lO^^'^eV^ theory 
and experiments agree reasonably well, with the exception of the Gallium experiments, in 
which case an agreement is obtained roughly within two standard deviations. The allowed 
values close to Ar?7,2i — 10~^°eV^ are in agreement with those obtained by two-generation 
analyses ||22|,|l9|. 



V. CONCLUSIONS 

Assuming maximally-mixed three-generations of neutrinos, we have analyzed the recent 
data on atmospheric neutrinos obtained by the Kamiokande Collaboration, which presented 
a corresponding analysis in the two- neutrino scheme. 

The study is motivated by the indications of possible large mixing angle solutions for 
both the atmospheric and solar neutrino problems. In particular, in the recent Kamiokande 
analysis of the atmospheric neutrino data, the best fit values of sin^(2'(9) are shown to 
be almost unity. Both the sub-GeV and multi-GeV data can be explained with Aml2 — 
2 X 10~^eV^. In our three-generation case the range of values of Am^j agrees with that of 
the two-generation analyses, but the values of sm'^{2{}af3) are 0.89 {a , P = e , fi , t) because 
the maximal mixing angles are 35° instead of 45°, which is the case in the two-generation 
schemes. We conclude therefore that, besides this difference, the addition of the third 
generation of neutrinos does not change the range of Am32 obtained in the two-generation 
analysis. 

Our solutions are based on vacuum oscillations. For the range of Am32 that we obtained, 
the medium effects of the Earth are negligible. Downward-going neutrinos do not travel 
enough in matter and for upward-going neutrinos the large neutrino mixing and the value 
of A?Ti32 make the MSW effect inoperative. 

We have also checked the solar neutrino problem in the same scheme. One definite 
prediction of the maximally mixed three-generation model is that the observed event rate 
for the Gallium experiments cannot be larger than 80 SNU. With exception of the Gallium 
case, the scheme provides a good fit of the data, assuming vacuum oscillations with values of 
Am^i — 10~^°eV^, similar in magnitude to those obtained in the two-generation analyses. 



Therefore, we conclude that the maximally-mixed three generations of neutrinos are 
capable of explaining both the current atmospheric and solar neutrino problems (though 
marginally in the case of the Gallium solar neutrino experiments). Improved data in the 
future will decide the fate of this scheme. 

Finally, the obvious question is: Is it possible for completely non-degenerate neutrinos 
(with the standard mass hierarchy mi <^ 1712 <^ WI3) to have maximal or, at least, very large 
mixing? The conventional wisdom tells us that mixing angles are determined by lepton mass 
ratios. In the quark sector, the Cabibbo angle is well reproduced by the formula tan^ ^^ = 
fTT'd/'^s- In the same spirit, the neutrino mixing angle, in the case of two generations, is 
expected to be determined by the ratios me/rn^ and m{i'i) / m{v2) ■ If wi(z/i) <^ rn{h'2), one 
expects the neutrino mixing angle to be very small. 

However, this conclusion is not valid when one invokes the see-saw mechanism. The pos- 
sibility that the see-saw mechanism may enhance lepton mixing up to maximal was discussed 
in Refs. PU| , P(]| . When the see-saw mechanism is invoked with more than one generation of 
neutrinos, the mixing angles for the resultant light Majorana neutrinos become dependent 
upon the masses of heavy right-handed Majorana neutrinos which are usually provided by 
GUTS models. In particular, when certain relationships among the mass parameters in the 
original Dirac and heavy right-handed Majorana mass matrices are satisfied, the mixing 
angles for the light Majorana neutrinos can substantially be enhanced, even in the case of 
the usual hierarchy m(z/i) ^ r7i(z/2) ^771(2/3). 

A full analysis of the atmospheric neutrino problem with arbitrary 3 mixing angles will 
be given elsewhere. 
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FIGURES 
FIG. 1. Separate fits of the sub-GeV (a) and multi-GeV (b) data for several different values of 

FIG. 2. Values of x^ for the separate fits of the sub-GeV (a) and multi-GeV (b) data and the 
combined fit (c). 



FIG. 3. Solar neutrino event rates as functions of Airily. The horizontal solid lines are the 
central values of the experiments with the dashed lines representing one standard deviation. The 
oscillating solid lines represent the theoretical values, with the dotted lines indicating one standard 
deviation resulting from the errors in the SSM. 
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FIG. 3. Solar neutrino event rates as functions of ^m\-^^. The horizontal solid lines are the 
central values of the experiments with the dashed lines representing one standard deviation. The 
oscillating solid lines represent the theoretical values, with the dotted lines indicating one standard 
deviation resulting from the errors in the SSM. 
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